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Current theories of the origin of the Universe, including string theory, predict the 
existence of a muhiverse containing many bubble universes. These bubble universes will 
generically collide, and collisions with ours produce cosmic wakes that enter our Hub- 
ble volume, appear as unusually symmetric disks in the cosmic microwave background 
(CMB) [1] and disturb large scale structure (LSS). There is preliminary observational 
evidence consistent with one or more of these disturbances on our sky [2]. However, 
other sources can produce similar features in the CMB temperature map and so ad- 
ditional signals are needed to verify their extra-universal origin. Here we find, for the 
first time, the detailed three-dimensional shape and CMB temperature and polariza- 
tion signals of the cosmic wake of a bubble collision in the early universe consistent 
with current observations [3]. The predicted polarization pattern has distinctive fea- 
tures that when correlated with the corresponding temperature pattern are a unique 
and striking signal of a bubble collision. These features represent the first verifiable 
prediction of the multiverse paradigm and might be detected by current experiments 
such as Planck [4] and future CMB polarization missions. A detection of a bubble 
collision would confirm the existence of the Multiverse, provide compelling evidence for 
the string theory landscape, and sharpen our picture of the Universe and its origins. 



1 Introduction 

The possibility of observing cosmic bubble collisions has recently received a considerable 
amount of attention (see, e.g., [5] for a recent review). Such collisions are a generic 
prediction of multiple- vacua models like the string theory landscape, and an observation 
of one would fundamentally alter our understanding of the cosmos at large scales. In 
these models the observable part of the Universe is contained within a bubble that 
formed as a result of a first-order phase transition from a parent false vacuum. Collisions 
with other such bubbles produces a special wave, that we term a cosmic wake, that 
propagates into our bubble and affects the spacetime region to the causal future of the 
collision (see Fig. 1). 

To a remarkable extent these dramatic events can be analyzed analytically and in 
a model-independent fashion. In this work we determine the precise pattern of temper- 
ature and polarization of Cosmic Microwave Background (CMB) photons induced by a 
cosmic wake using the full solution to the cosmological Einstein-Boltzmann equations. 
The results, while consistent with a previous analytic approximation [6], are remark- 
able: collisions can produce a unique and highly characteristic polarization signal — a 
"double peak" in the magnitude of the polarization as a function of angle. This dou- 
ble peak, and more generally the large-scale azimuthally symmetric polarization and 
temperature pattern produced by cosmic bubble collisions, serves as a true smoking 
gun for their detection. Moreover, in an interesting regime of parameters this vital 
corroborating signal in CMB polarization can be as easy or potentially even easier to 
detect than the signal in CMB temperature. 

Our results also determine the evolution of density and velocity perturbations after 
decoupling, which opens the door for future work quantifying the effect of the cosmic 
wake on large scale structures. 

2 Background 

While the detailed physics of cosmic bubble collisions depends on high-energy micro- 
physics, to a surprising extent the physics relevant to the late-time cosmic wake that 
affects observables such as the CMB radiation and Large Scale Structure (LSS) appears 
to depend on only a few simple parameters. This occurs because the collision between 
two bubbles preserves much of the symmetry inherited from bubble nucleation events, 
and the inflationary period in our Universe subsequent to any collision rapidly erases 
all but a few characteristic features. 

In a broad class of microphysical models there are only four parameters relevant to 
the effects of a collision on the observable Universe at the end of inflation. Of these, one 
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Figure 1. Spacetime conformal diagram of a cosmic bubble collision. Each point represents a 
two-dimensional hyperbolic slice of space with a radius of curvature that depends on position 
within the diagram. The physics in regions A, B, and C is determined by the conditions on 
the surface labelled "initial cond." inside the cosmic wake and by the conditions outside the 
cosmic wake. To a good approximation, the physics in regions B and C depends only on the 
left-moving part of the inflaton perturbation. 

probes the underlying microphysics, while the others reflect the initial conditions (the 
position and time of the nucleation of the colliding bubble relative to ours). These four 
parameters determine the direction, distance to, and strength of the collision-induced 
cosmic wake in our Hubble volume, and therefore the location, size, and intensity of 
the affected disc in the CMB sky. 

Because the primordial plasma supports acoustic waves, the cosmic wake can prop- 
agate into our Hubble volume. The most dramatic effect is the appearance of two sharp 
peaks in the magnitude of E-mode polarization as a function of radius: there is a double 
concentric ring of strongly polarized light outlining the affected disk in the CMB sky. 
These peaks are a position-space manifestation of the acoustic physics responsible for 
the more familiar peaks observed in the angular temperature power spectrum of the 
CMB. 

2.1 Scalar Perturbations after Inflation 

We first calculate the expected curvature perturbations resulting from cosmic bubble 
collisions that persist after a period of slow-roll inflation. These perturbations subse- 
quently evolve into the cosmic wake. 

The spacetime metric and scalar-field configuration describing the collision of two 
Coleman-de Luccia thin-wall bubbles expanding in a parent false vacuum preserves an 
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5*0(2, 1) group of isometries [7-10]. In suitable coordinates the collision takes place at 
one instant of time everywhere along a hyperbolic surface EI2 in space, and the effects 
of the collision causally influence the spacetime region to the future of the surface at 
that instant. 

During inflation the background metric describes, up to slow-roll corrections, an 
approximate de Sitter spacetime with Hubble constant Hi and may be written in the 
form ds^ = -(1 + Hft^y^dt^ + (1 + Hjt'^)dx^ + t'^dMl where dM^ = dr^ + sinh^ rdip'^ is 
the metric on the surface EI2. Since the collision preserves the 5*0(2, 1) symmetry its 
effects are uniform on the IHI2. 

However, after N e-folds of inflation the coordinate Hjt oc becomes expo- 
nentially big and the hyperbolic curvature small. Current constraints indicate that 
-Ro ^ Wc/Hq, or equivalently a/H^ = {RqHq/c)^^ < 0.1, where Rq is the radius of 
curvature of the Universe today and Qk is the effective curvature density [3], so it is 
a good approximation to ignore the negative spatial curvature and treat the hyper- 
bolic surface EI2 as planar (which we do for the remainder of this work). We also 
approximate the inflationary spacetime as de Sitter space. Curvature and higher-order 
slow-roll effects induce small corrections to the result calculated here. Finally, we use 
the thin-wall approximation for the collision bubble — that is, we treat the region af- 
fected by the collision as having a sharp boundary. Taking into account finite thickness 
will have the effect of smearing out very small scale features in the signal, and sufficient 
smearing would replace the double peak with a single peak. The extent to which this 
happens is a model-dependent question and could provide further details about the 
microphysics — specifically, if a distinct double peak is observed it indicates that the 
ratio of the Hubble scale during infiation Hj to the Hubble scale of the parent false 
vacuum Hp satisfies {Hi/HpY < 10"^ [11]. 

In the planar approximation and including linear perturbations in conformal New- 
tonian gauge, the metric takes the form 

ds^ = (Hit)-^ [-(1 + 2^)dT^ + (1 - 2^){dx^ + dy^ + dz^)] , (2.1) 

where Hj is the Hubble parameter during infiation, r < is the conformal time, the 
spatial coordinates are chosen so that the solar system's comoving location is a?' = 0, 
and $ and \1/ are the Newtonian potentials (which are equal during infiation). In 
the approximation that the collision surface is planar, the region of spacetime affected 
by the collision satisfies x > — r + Xc for some constant Xc that depends on when 
the collision occurred. We focus on collisions with a region of causal infiuence that 
intersects the visible part of the surface of last scattering because these have the best 
potential for detection [1, 6, 12]; since the Earth is at x = and Hj\t\ ^ 1 at the end 
of inflation, we take Hj\xc\ <^ 1. 
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Given that slow-roll inflation took place, pre-infiationary inhomogeneities became 
small at some time early in inflation. From that time on we may evolve perturbations 
using standard linear cosmological perturbation theory [13]. To lowest order in the 
slow-roll parameters, gauge-invariant perturbations to the inflaton fleld (p = (po + 
obey the equation of motion of a free scalar in de Sitter space [13] 

dl5(t) + 2HdM-W^5(t) + 0{t,r])=Q, (2.2) 

where 'H(r) = — 1/r + (^(e,//) is the conformal Hubble constant, and e,?7 are the 
standard slow-roll parameters. We neglect terms of 0{e,ri) and work at zeroth order 
in the slow roll expansion. 

Because the collision perturbation is constant on the nearly-planar surface the 
general solution to (2.2) can be written in closed form [1] 

50(r, x) = g{T + x)- T~g'{T + a;) + fir - x) - rf'{T - x), (2.3) 

where / and g are arbitrary functions of one variable (and /' and g' their derivatives 
with respect to their argument). To proceed further we need the initial perturbation 
6(f) produced by the bubble collision event. Without a speciflc model for the underlying 
microphysics this cannot be computed but, to a good approximation, the post-inflation 
predictions of any such model can be characterized by four parameters. 

Consider the region near the edge of the collision lightcone x = —ti + Xc at some 
time Ti ~ near the beginning of inflation when linear perturbation theory is 

valid (labelled "initial cond." in Fig. 1). Since the perturbation 50 is non-zero only for 
X > —Ti + Xc, from (2.3) we must have g{T + x — Xc) = gij + x — Xc)0(r + x — Xc) and 
fij—x) = /(— r+x— Xc+2rj)6(— r+x— Xc+2rj), where 9 is the Heaviside step function 
and /, g are functions of one variable. The / terms in (2.3) are rightmoving excitations. 
By the end of inflation Hj\Tf,\ ^ 1, they are non-zero only for x — Xc > — 2rj ~ 2HY^ 
but, as explained above, we are only interested in \x — Xc\ ^ Hy^ and these terms are 
not relevant for cosmological observables [c.f. Fig. 1, speciflcally regions B and C). 

At the end of inflation, the inflaton perturbation is [11] 

oo 

6(j){T,,x) ^ g{x - Xc)e{x - Xc) =M^{a„i77(x-Xe)"}e(x-Xc), (2.4) 

n=0 

where we have dropped the term in (2.3) proportional to t^, ^ Hy^, and expanded g 
in terms of dimensionless coefficients a„ and a constant M with dimensions of mass. 

The model-dependent initial conditions (50(rj,x) and Scjyiji^x) determine the coef- 
ficients an- As is clear from (2.3), any regular perturbation 50 gives ao = 0. For an 
0{1) perturbation at an early time Xj ~ Hj^ one generically expects the dimensionless 
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coefficients a„ < 1 forn > 0. But as we are only interested in the region Hi\x—Xc\ ^ 1, 
absent fine-tuning of tlie tlie fiiglier terms in (2.4) will be negligible. Thus, by the 
end of inflation, the leading perturbation reduces to just the linear term [1, 11] 

6(f) = Mai{x - Xc)Q{x - Xc). (2.5) 

This expression contains two parameters {Mai and Xc)- We have chosen coordinates 
so that the x direction is the axis of symmetry of the collision; i.e. there are two angles 
parametrizing the direction towards center of the collision. As mentioned above, these 
four parameters fully characterize the effects on the CMB of a generic collision. 

2.2 Curvature Perturbations after Inflation 

The infiaton perturbation 6(f) determines the conformal Newtonian gauge gravitational 
potential perturbation $. Using the Einstein and scalar equations we can show (2.5) 
leads to a potential perturbation 

v v 

$ = -y^(^ + x)e{T + x)^ -y Maix0(x), (2.6) 

where V = dV{(f>)/d(f> is the slope of the inflationary potential. 

The initial conditions for cosmology after inflation are best expressed in terms 
of the comoving curvature perturbation ( that is conserved on super-Horizon scales 
[13, 14]. We find 

C = 7T 7^ + ^ = 7, T— + T + x]e{r + x)^ AxG X , 2.7 

3 1 + w 3 \l + w / 

where ti? ~ — 1 is the equation-of-state during inflation and A = —aiMV'{2 + w)/{3 + 
3w) ~ aiMV/V is a constant that sets the amplitude of late-time effects. The final, 
approximate equality is valid late in the inflationary epoch. 



3 Results 

Using the superhorizon curvature perturbation Q near the end of inflation and a set 
of linear transfer functions Ax the perturbed distribution of an observable X at later 
times is 

Ax(x,n,r) = j d^k e^'^'^Axik = \k\,n,r)Q{k), (3.1) 

where Ax(x, n, r) is the local value of an observable X (e.g. 6T/T for the photon 
distribution) at position x, and time r in the direction n on the sky, and k is the 
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comoving wave-vector. Angular moments of Ax(x, n, r) determine, for example, the 
density and velocity as a function of position, while Ax(0, n, r) is the angular anisotropy 
at the location of the Earth. 

The transfer functions encaspulate the full evolution of the coupled multi-component 
fluid and gravitational system and to linear order depend only on the background 
cosmology — hence they are independent of the collision perturbation. For the cosmo- 
logical background we choose best-fit values from the 7-year data release of WMAP 
[3]. We compute the transfer functions using a customized version of CAMB [15], 
and our own code to numerically evaluate the Fourier transform in (3.1) and extract 
observables. 

3.1 Temperature 

In the Sachs- Wolfe approximation [16], the CMB temperature anisotropy is determined 
by the curvature perturbation at decoupling: 6T{9,(t>)/T oc ({x,Tdc), where = D^^ 
is the distance to the last scattering surface and the angular coordinates are chosen 
so that X = DdcCosO. In this approximation the collision perturbation (2.7) gives rise 
to a very simple temperature perturbation: ST/T oc Q{9c — 0) {cos 9 — cos 6'c), where 
6c = cos~^ Xc/Ddc is the angular radius of the spot [1, 6]. 

However, the perturbation (2.7) is constant only on super-horizon scales after the 
end of inflation. To take into account the full evolution between reheating and decou- 
pling, we instead use (3.1). The results are presented in Fig. 2, with = pointing 
along the x-axis (towards the collision) and \5T /T\ normalized to 10~^ at the center of 
the spot. We present a range of Xc on the figure: the ?/-axis is the approximate angular 
radius of the spot O^. ^ cos~^ xd Ddci and the x-axis is the angle from the spot center. 

For all Xci we find the temperature perturbation has the largest magnitude at the 
center of the spot ^ = 0, decreases linearly in cos^, and then smoothly transitions to 
a constant at 9 = 9c (with no edge or discontinuity). The width in 9 of the transition 
region is roughly equal to the angular separation between the two peaks in polarization 
discussed below. These results are fully consistent with the approximate analytic results 
of [1, 6]. Whether the spot is hot or cold encodes microphysics of the inflaton and the 
collision bubble, see [1]. 

3.2 Polarization 

The CMB is linearly polarized due to Thomson scattering of CMB photons off free 
electrons [17, 18]. This scattering occurs primarily at redshifts around decoupling 
{zdc ~ 1100) and reionization {zre ~ 10). Scattering produces a net polarization 
proportional to the quadrupole temperature anisotropy of the radiation incident on 
the electron. Scattering at decoupling depends on well understood atomic physics [19], 
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Figure 3. Polarization for (a): a spot with angular radius 9c ~ 10°, and (b): a 

spot with angular radius Oc ~ 60°. Note the double peak in (b), with the indicated region 
magnified in (c). The spots are normalized so that \5T /T\ = 10~^ at the center. 

but there is significant uncertainty in the reionization history of the Universe. However, 
our key results are not very sensitive to different reionization models [6], and we choose 
a simple model of single reionization as our fiducial case [3]. 

The polarization of a transverse electromagnetic wave with intensity tensor Jjj can 
be characterized by the Stokes parameters Q and U: 

Q = {hi - /22)/4, U = /12/2. (3.2) 

Choosing coordinates so that the colhsion spot is centered on 6* = guarantees that 
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Figure 4. Density plot of the polarization for spot radii from 9c = 5° to 90°. The plot 

is normalized such that each spot has a temperature contrast \5T/T\ = 10~^ in the center. 

the Stokes parameter f/ = {i.e. the polarization is purely ii^-mode, as expected for 
a scalar perturbation). The polarization pattern is radial (azimuthal) for a cold (hot) 
spot with Q > (<) 0. For a discussion of CMB polarization in general see [20]. 

Thanks to the azimuthal symmetry, we can parametrize the CMB polarization due 
to a single collision hy Q = Q{9) alone. However Q{9) differs qualitatively for small 
{Oc < 12°) and large spots. We display two cases in Fig. 3; 9c ~ 10° and 9c ~ 60°. 
The large spot has a distinctive "double peak" structure; two rings of sharply enhanced 
polarization concentric with the edge 9 = 9c oi the affected disk in the temperature 
map. This structure is characteristic of bubble collisions and appears for any spot with 
9c > 12°. 

In Fig. 4 we display the results for all collisions with 5° < < 90°. The x-axis is 
9c, the y-axis is |(5(6')|. 

3.2.1 Why a double peak? 

The double peak in Fig. 3 and Fig. 4 is a very striking feature, but the physics that gives 
rise to it is simple. At the end of inflation, the perturbation (2.6) in the Newtonian 
potential is $ ~ xQ{x). Between the end of inflation and decoupling, a linear- in- 
X perturbation (such as $ = Cx) remains so, because the evolution equations for 
cosmo logical perturbations are all second order in spatial derivatives. By contrast, the 
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Figure 5. The last scattering surface and the cosmic wake, showing the region of the CMB 
temperature sky affected by the collision and illustrating the geometric origin of the double 
peak separation angle. 

"kink" at a; = spreads at the sound-speed Cg < c/a/3 in the plasma, and evolves 
into a smooth profile confined within the sound horizon of a; = until matter-radiation 
equality. After equality the speed of sound is effectively zero and there is no further 
spreading. However at the edge of the affected region at any time, the first derivative 
of the perturbation is still close to discontinuous. 

An electron with a last scattering sphere that intersects one of the edges of this 
region will scatter a relatively large amount of polarized light, because the quadrupole 
moment is sensitive to the second derivative of the temperature distribution. The 
two edges of the region at decoupling should be separated by twice the sound horizon 
2rsh ~ 306 Mpc, or Tsh/D^c. ~ .0109 [3] (see Fig. 5). Therefore the angular separation 
between the two edges should be 6s ~ cos~^ (xc — Tsh) / Ddc — cos~^ (xc + Tsh) /Ddc, 
which agrees well with our numerical results. For 9c ~ 90° this width is 9s ~ '29 sh = 1-3°. 
For smaller disks, 9s is larger, and for 9c < 12° the inner peak disappears entirely. 

3.3 Density 

After matter-radiation equality, the sound speed of perturbations is approximately 
zero and the edge of the cosmic wake should remain at fixed x. Due to the growth of 
perturbations during matter domination, the amplitude of the perturbation at the edge 
of the wake today should be amplified by roughly a factor of Zdc ~ 1100 from its value 
at decoupling. Since a Newtonian potential of the form $ = Cx9{x) corresponds to a 
5-function density perturbation at x = 0, one expects a planar sheet of over- or under- 
density (for hot or cold spots in the CMB, respectively) centered at x ~ DdcCos9c, 
with a thickness set roughly by the sound horizon Vgh at equality. 
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X (Mpc) X (Mpc) 

Figure 6. Evolution of the cold dark matter density contrast 5cdm in conformal Newtonian 
gauge, from redshift z = 10, 000 to the present. The left panel shows the continuous evolution, 
while the right panel shows two snapshots at z = 10,000 (top) and z = (bottom). The 
cosmic wake is normalized so that 5T/T = — 10""^ (a cold spot) at the center of an affected 
disk in the CMB with angular radius 6c ~ 80°. 

These expectations are confirmed numerically in Fig. 6, where we display the cold 
dark matter density contrast in conformal Newtonian gauge (the baryon density looks 
similar) . The presence of such a mass sheet could lead to signatures in large scale struc- 
ture or lensing surveys which — combined with CMB data — could further corroborate 
the discovery of a bubble collision. We leave this investigation to future work [21]. 

4 Detectability 

In this section we estimate the degree of detectability of the polarization signal for spots 
of varying size and brightness and for a selection of current and future experiments. 
Both the temperature and polarization signals are circularly symmetric. We orient our 
coordinates such that the center of the spot is at the pole 6* = 0; a coordinate rotation 
can be used to center the spot anywhere else. 

To get an estimate of the detectability threshold we will fit to a one parameter 
model (which we label A) for the amplitude of the temperature spot at its center. 
We analyze the fit for temperature (T) alone, E-mode polarization (E) alone, and 
E-mode polarization combined with information from the temperature map and cross- 
correlation. We work in Z-space where the noise covariance matrix is diagonal. The 
circular symmetry allows us to integrate over the azimuthal angle (p, which in /-space 
means the multipoles do not contribute for m j^O. 
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For each estimate we generate a full sky of CMB fluctuations in T and E con- 
sistent with WMAP-7 concordance parameters [3] as well as a noise realization for 
each detector up to Imax = 2000. We then add a collision spot to this and compute 
the relative likelihood function (which amounts to a simple this one parameter 

model) to determine if the spot is detectable. We normalize such that for a given 
spot a null detection is A = (i.e. consistent with Gaussian fluctuations) and a per- 
fect detection is A = 1. As a sample we look at the projected sensitivities for three 
experiments: a current satellite mission (Planck [4]), a planned balloon mission (SPI- 
DER [22]) and a future satellite mission (CMBPol (EPIC-2m) [23]). We use the same 
predicted sensitivies, beam- widths, weight per solid angle (w), and observing time as 
[24]. For brevity we analyze just these three, but any experiment which measures the 
temperature and/or polarization over a significant fraction of the sky — such as CLASS, 
EBEX [25], PIXIE [26], ACTPol [27], and SPTpol [28]— has the potential to detect a 
cosmic wake. 

Since we are fitting for the aim themselves rather than the power spectra Q, a cut 
sky introduces straightforward but bothersome complications (the spherical harmonics 
are no longer orthonormal) . To get a simple idealized estimate we use a full sky for 
each. For the satellite experiments this amounts to ignoring foreground effects in our 
estimates. For partial sky-coverage balloon missions we are assuming that the Gaussian 
fluctuations will be similar in the region of the sky not measured. A full treatment of 
these effects would likely decrease the quoted sensitivities by a modest factor. 

For each experiment we do a combined analysis across all bands. Table 1 sum- 
marizes these numbers (we have assumed the sensitivity to Stokes / is \/2 that of Q 
and U). For each spot radius 6c and experiment we find the minimum \6T/T\ at the 
center of the spot necessary such that a detection A = 1 is 3a away from zero. This 
represents a conservative estimate of the detectability threshold in temperature and 
E-mode polarization. We summarize our findings for temperature and polarization in 
Fig. 7. For polarization we can also choose to use information from the cross-correlation, 
which increases the sensitivity by allowing us to eliminate the correlated part of the 
contribution from Gaussian fluctuations. We compute a range of sensitivities using as 
much information from the cross-correlation as possible to none which is displayed in 
Fig. 7(b) by the filled region for each experiment. 

We can make some general observations. For small spots {6c < 10°), a spot with 
\6T/T\ > 4 X 10~^ is likely observable in temperature, but would have to be brighter 
than ~ 2 X 10~^ to be observable in polarization. A preliminary search for bubble 
collisions in the temperature map has been recently carried out [2, 29] using the WMAP- 
7 dataset. While the analysis detected several anomalous features, it rules out collisions 
with \6T/T\ > 10~^. Thus, while the potential still exists to observe small spots in 
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Table 1. Projected detector parameters 




Figure 7. 3a detectability thresholds for (a) temperature and (b) E-mode polarization. In 
each, Planck is red, SPIDER is blue and CMBPol is green. The thickness of the lines in (b) 
corresponds to using information from the cross-correlation to improve the measurement. 

the temperature map, correlating them with polarization signals is likely outside the 
capabilities of current or near future experiments, at least at the 3a level. If we lower 
our threshold for correlation it may still be possible to obtain some evidence for a 
collision in the polarization map, though not likely significant enough to exclude other 
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potential explanations. 

The situation is considerably more promising for larger spots. With 6c > 25° much 
fainter spots are detectable in polarization, and very large spots could be detected 
almost as easily in polarization as temperature. The reason for this is clear from 
Fig. 3. Larger, fainter spots have stronger Q (or E) -mode polarization over a larger 
angular range and thus more easily can be seen over the Gaussian fluctuations and 
detector noise. In addition, spots larger than approximately 12° have the distinctive 
double peak feature on smaller scales. 

This angular size dependence is particularly important because the size distribution 
of collision spots — dN{Oc) oc sin^c'^^'c — is a robust prediction of the theory [12]. The 
analysis of [2, 29] was restricted to spots with Oc < 11°, which, according to this 
distribution, are much less common than larger spots. 

5 Conclusions 

We have solved numerically for the real-space evolution of the cosmic wakes produced 
by collisions between bubble universes, and determined the distinctive temperature and 
polarization patterns these cosmic wakes imprint upon the CMB. These patterns have 
a circular symmetry that reflects the near-planar symmetry of cosmic wakes, and the 
polarization pattern can have a distinctive double-peak structure arising from propa- 
gation of acoustic waves in the primordial plasma. We have estimated the detectability 
of cosmic wakes for several current and future CMB experiments in both temperature 
and polarization. Holding fixed the temperature anomaly at the center, increasing the 
radius of the spot makes it easier to detect with polarization, but harder using temper- 
ature. The detection of a cosmic wake would show our observable Universe is one part 
of a larger multiverse, support the idea of the string theory landscape, and constitute 
a groundbreaking discovery of the nature of the big bang. 
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